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ABSTRACT 

We report HEAO-1 A- 2 high time resolution X-ray observations of the 
X-ray emitting Seyfert I Galaxy NGC 6814. In sharp distinction with a sample 
of over 30 active galactic nuclei (Tennant et al. 1981 > this object showed 
strong X-ray variability on timescales less than 3 hours. The mean flux on a 

timescale of 90 minutes varied by a factor of - 2.5 corresponding to AL 'lx 
43 

10 ergs/sec. An autocorrelation analysis shows a characteristic time for 
+60 

variability of 100_ 2 ^ seconds. There is no indication of spectral variability 
with an upper limit on a change in the power law spectral index 
of | AT | < .37, for a factor two change in intensity. We consider in detail 
the constraints of such rapid variability on a wide variety of X-ray source 
mechanisms . 

Subject headings: galaxies: individual - galaxies: Seyfert - 

X-rays: sources 
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X. INTRODUCTION 

X-ray observations of the tine variability and spectra of active 
galactic nuclei (AON) relate directly* in the present conventional wisdom, to 
the engine responsible for the enormous luminosities of these objects, Xn 
particular, observations of rapid X-ray variability plaoe strong constraints 
on the sise of the "central engine", information which is difficult to obtain 
from observations in the radio, infrared, optical or ultra-violet bands. 

Various models of the energy generation mechanism in AGN's have centered 

(Fabian 1979) on the existence of a compact object in the nucleus of these 

galaxies. Indirect arguments and various scaling laws indicate that this 

putative central object has an effective scale sise of less than 10^ cm 

(Fabian 1979). Such an object vould be revea'ed by characteristic time 

variability of less than one day. A general constraint on the fastest 

timescales ejected has been proposed by Fabian and Rees (1979), building on 

work by Cavallo and Rees (1978). If a spherically symmetric, homogeneous, 

cloud of stationary matter is converted into radiation with efficiency n, the 

shortest observed timescale for the conversion occurs when the optical depth 

of the cloud is unity. This gives a minimum timescale At . ' 5 L . ,/n sec, 

where l 43 is the luminosity in units of 10 43 ergs/sec. Typical AGN's have 

^43 — 1 and, if n is as large as 0.1, timescales as short as SO sec are 

possible. On the other hand models of accretion disks (Shapiro, Lightman and 

Bardley 1976) indicate that most of the X-ray luminosity would come from a 

GM 

region of sise R ){ »■ 10 —j . If the luminosity of these objects is Eddington 

c 

limited then R^ « 10’ 1 l^ 3 cm for a characteristic time of * 4 L 43 sec. 

Clearly X-ray observations capable of measuring timescales in the 1-10,000 sec 
time range will sample the predicted characteristic times of the "central 
engine" in AGN's. Previous X-ray observations of time variability in AGN's 
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have been prevented from sampling such short timescales by a combination of 
law source flux, low effective area and relatively high internal background* 
For only a few objects, NGC 4151 (Tananbaum et al. 1978), Cen-A (Delvaille, 
Epstein and Schnopper 1978), and NGC 4051 (Marshall et al. 1981) are there 
reported timescales less than 6 hours* In addition, Tananbaum (1980) has 
reported timescale of ■* 6 hour for NGC 6814. 

In this paper we report the detection of short timescale fluctuations 
(t < 1 hr) in the X-ray flux from the Seyfert I galaxy NGC 6814* This 
observation is the result of a larger program to look for short term 
variability from AGN's, the results of which will be reported elsewhere 
(Tennant et al. 1981). 

II. OBSERVATIONS AND DATA ANALYSIS 

The GSFC detectors on the HEAO- 1 + satellite have been described in 

+ The A2 experiment on HEAO- 1 is a collaborative effort led by E. Boldt of GSFC 
and G. Garmire of CIT, with collaborators at GSFC, CIT, JPL and UCB. 

detail by Rothschild et al. 1979. Because of the peculiarity of our results 
and the importance of ensuring their validity we describe in Appendix A the 
analysis procedure used to derive the NGC 6814 data. 

Briefly, the observations of NGC 6814 were performed by the HEAO-1 
satellite experiment A-2 on April 28 of 1978 from 0 hours UT to 8 hours UT 
while the satellite was in the pointed mode. Counting rates were returned 
every 5.12 sec from 8 "independent" counters. Background was determined from 
deep exposures of "source free" regions of the sky where no sources of flux > 

2 x 10 ergs/cnr/sec were seen. Thr net counting rate was corrected for the 
effective exposure of the detectors to the source. The rates from the 
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different detectors could them be sunned together to construct e "light 
curve", with the best available statistics • 

In Appendix A we discuss the wide variety of checks made in order to 
ascertain that the strong variability seen in HOC 6814 was not a result of 
instrumental problems* We conclude that NGC 6814 is the source of the 
variable X-ray flux detected with HEAO-1 A-2. 

III. RESULTS 

Examination of Fi ur 1 shows that the flux from NGC 6814 is strongly 

variable. On an orbit-b* -orbit (1.5 hour) basis the flux varies by a factor 

2 

of * 2 .,5 (Table 1). In Table 1 we also show the x calculated on the 

assumption that the source was constant during each individual orbit. For all 

2 

but one of these orbits the value of x allows us to reject the hypothesis of 
source constancy on timescales less than 1.5 hours at > 99.99% confidence. 

This is clearly shown in Figure 2 xn which a prominent flux change on a 
timescale of less than 500 seconds is evident. 

Having established source variability on timescales less than 1.5 hours 
we now attempt to characterise the nature and timescale of this variability. 

A. Timescale of variability 

The autocorrelation function (ACF) was confuted for the first and second 
orbits using a 10.24 sec bin size. In Appendix B we discuss how the ACF was 
generated and how errors in the reported timescales were determined. The 
inferred characteristic timescale of the source intensity and the errors are 
given in Table 2 for these orbits and the entire data stream. All the data 
are consistent with a characteristic time of ' 100 seconds. There was 
insufficient data for orbits 3-5 for a meaningful autocorrelation function to 
be calculated separately. In Figure 3a we show the ACF for the complete data 
set (orbits 1-5 inclusive). The characteristic time for variability is the 
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slope of the natural logarithm of the ACF (Figure 3b). In panels 3c end 3d 
are shown, respectively, the "detrended" ACF function and its natural 
logarithm for the sms data (see Appendix B for details). The timescale of 
130 seconds seen when the data are "detrended" confirms that this is a 
characteristic time for the entire data set. 

An exponential ACF could be interpreted in terms of a shot noise 
model. However, we note that many other processes give similar ACF' 8 
especially when there is a large amount of statistical noise in the dattk. In 
addition, it is quite clear from Figure 1 that the shot parameters change from 
orbit to orbit. A final note is that a shot noise model is only a method of 
parameterizing the data and does not necessarily have any physical 
significance. 

B. Limits on Spectral Variability 

He have searched for spoctral variability in four different ways. 1) 

He have computed the mean spectrum on an orbit by orbit basis. 2) He have 
sorted the data on the value of their flux and have computed a mean "hardness" 
ratio for each value of the flux. 3) He have computed the cross-correlation 
function between two energy bands. 4) He have looked for spectral 
variability during a given event. Bach of these methods is most sensitive to 
different types of spectral change as discussed below. 

1. Orbit-by-Orbit Variability 

The best fit spectra on an orbit-by-orbit basis assuming a power law 

model of photon index T and equivalent column density of hydrogen in the line 

of sight N|f are shown in Table 3 (all errors are 1o). These values are to be 

compared with the flux weighted average for the entire observation of T * 

1.73±.22, N|| ” 4«3+2*o x *° 22 cn ~ 2 (Mushotzky et al. 1980) (90% confidence 

2 

error) . The value of % computed on the assumption that T and n h for each 





i 


I 

i 
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orbit is ths same as the flux weighted mean is 6*8 for 10 parameters of 
interest* This indicates that, despite a factor of 2 change in intensity 'n a 
timescale of 1.5 hours, the time averaged spectrum did not vary significantly. 

If we assume that N H did not vary during this observation (consistent 

with che values in Table 3) and fix N H at the value derived for the entire 

observation, the formal error in r decreases. In Table 4 we show the best fit 

values of T for each orbit under this assumption. There is a slight 

indication that the source spectrum is marginally flatter when the flux is 

hitter. In Figure 4 we plot T vs the 2-20 kev flux on an orbit-by-orbit basis 

and find that a linear relation of the form T • mF+b (where F is the flux) 

+ •022 

gives an acceptable fit with m - -.054 .*032 * I * 1 * 8 trend is significant at 

only 2.5o however. Alternatively we place a 30 upper bound of AT < .37 on an 
orbit-by-orbit basis. 

2. Limits on Spectral-Intensity Correlation 

We have accumulated data at B > 6 keV ( and 5 < 6 keV (E Lqw ) 
on an 82 second basis. We then sorted the data into 10 intensity bins on the 
basis of E Low - In Figure 5 we plot the average of E High whenever E Lqw was 
within the prechosen bin limits, along with lines which correspond to .power 
law spectra of slopes T ■ 2.0, 1.75 and 1.50 with N H fixed at 4.3 x 10 22 
cm* 2 . The data are all consistent with the r * 1.75 line. We note that each 
point in Figure 5 represents an independent measure of the power law index. 

We conclude that there is no intensity related spectral variability with an 
upper limit on a change in spectral index of LT < .30. (This method of 
analysis could hide possible variability if the same intensity state has 
different spectral indices and if these values scatter about a given mean 
value which is also characteristic of the total time averaged spectrum. ) 

3 • Spectral Cross-Correlations 
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A defect of the prior two methods is the relatively large bin sises 
(30 minutes and 82 seconds) , necessitated by the relatively low flux* compared 
to the characteristic variability time of * 100 seconds* The 
cross-correlation function allows one to test the entire data stream on 
shorter timescales to see whether the softer photons systematically lead or 
lag the hard ones. We report here the result of two cross-correlation 
analyses* one designed to minimise the detectable timescale* and the other to 
maximise the difference between the selected energy bands* 

Cross correlation of the total flux from the NED (which is sensitive to 
photons with a mean energy of B * 5.46 keV for the NGC 6814 spectrum) with the 
total flux from HED3 (B - 8.67 keV) shows that the MED leads HED3 by 6.3 ± 3.3 
sec for a bin sise of 15.3 sec (see Appendix B). we also cross correlated 
photons of energy < 6 keV (E » 5.17 keV) with photons of energy > 6 
keV (B - 12.50 keV). For this analysis we used a bin sise of 41 sec since 
this made the errors due to Poisson noise similar to those in the MED vs HED3 
analysis. We found that the > 6 keV photons lead the softer photons by 4.1 ± 
9.6 sec. 

4. Spectral Variability of a Resolved Event 

The flux change shown in Figure 2 has sufficient statistics to 
examine spectral variability within a given event. As can be seen in the 
middle panel the effective spectral index T did not change (dr < .30) during 
this event despite the five fold change in source flux. These data also rule 
out variation In X-ray column as the origin of the short term variability. 

Due to the low flux the hardness ratios are binned in 82 second bins and thus 
we cannot, comment on spectral variability on the shorter timescale indicated 
in the previous section. We conclude that, for this one event* spectral 
variability does not occur on the characteristic time of intensity 
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variability* 

V* DISCUSSION 

Tha rapid variability obaarvad in NGC 6814 can place strong constraints 
on tha physical processes producing tha X-ray flux and tha environment in 
which it is produced. He shall discuss with respect to tha NGC 6814 data 
general constraints* and limits we can place on "standard" X-ray emission 
mechanisms. 

A. General Considerations 

The average luminosity of NGC 6814 is 1.1 x 10** h" 2 ergs/swc in the 

2-20 keV band where h is Hubble's constant measured in units of 50 

km/sec/Mpc. Since the amplitude of the variability is quite large we will 

assume that the change in luminosity is greater than 10 43 ergs/sec. If the 

X-ray spectrum extends out to 500 keV* as does the X-ray emission from Cen-A 

43 -2 

(Baity et al. 1981), then the average luminosity is « 4.7 x 10 h 
erg/sec . 

If NGC 6814 is powered by accretion onto a single compact object* then 
we can set a firm upper limit on the mass of this object by requiring that the 
fluctuation time be longer than the light travel time across the Schwar*. ''child 
radius . This gives 


M < 


c 3 At 

2G 


- 1 x 10 M 


for At ■ 100 sec. On the other hand* we find that 


M > 9 x 10 4 M„/h 2 


by requiring that the luminosity not exceed the Bddington limit. Without a 
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detailed modal we cannot restrict the mass range any further. As we will see 

though# a 10® M Q object is consistent with thermal Compton models. 

Using the Fabian and Rees (1979) relation given in the introduction 

with AL > 10 43 erg/sec and At . ' 100 sec# we find that n > 5%/h 2 . However# 

min 

if the variations extend out to 500 keV then n » 25%/h 2 . Xf one requires n < 
10% then h * 2. Conversely# if h * 1 then one must consider models that 
allow high efficiencies (Thorne 1974) or argue that one of the assumptions 
used in deriving the Fabian and Rees relation is violated. This could imply 
either non- spherical geometry# relativistic bulk motions# or that the source 
of the X-ray flux is not directly due to accretion. 

B. Emission Mechanisms 

1. synchrotron Radiation and Synchrotron Self-Compton Emission 

Xf the X-ray flux is produced by relativistic particles in a "steady 
state" type process via synchrotron emission# or synchrotron self -Compton 
(SSC) emission, then one requires that the particles "live" long enough to 
fill the emitting region. 

The lifetime of a relativistic electron against Cotton losses is 

T «* MS ' W here the energy of the electron is ym c 2 and u is the energy 

c yu 6 

density of the photon field. For L * 1.1 x 10 43 ergs/sec and a sise 
R * cAx ' 100 light seconds#u * 2 x 10 7 erg/cm 3 . 

Therefore t • y seconds. In order to produce X-rays via the synchrotron 
c 

2 ii 

process one needs y B > 6 x 10 (where B is the magnetic field in gauss). 

For "reasonable” values of B (< 10* gauss) one needs Y > 10* • Therefore for 

these particles t < 1 x 10~ 4 seconds. Xf the X-rays are produced by a single 
c 

Coup ton scattering off some photon population characterised by a 

2 

frequency then v x % Y (Jones# O'Dell and stein 1974). Xf we 
assign in the range 10 9 - 10 ,S Hs (radio to optical photons) #Y % 10 4 * 5 - 
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10 1 * 5 which gives Compton loss lifecimea of 10" 3 - 10"®*® seconds . 

All those timoscolos are much loss then ths light travel tins across the 
omitting region. Wo conclude that synchrotron and 88C models, without 
continual injection and/or reaccelereation and/or relativistic bulk motion* 
cannot work. 

Cavaliers and Morrison (1990) have considered a nodal in which the 
particles are continually roaccelerated and in which the radiation rate is 
controlled by the acceleration rate. Their model* for the parameters measured 

A 

in this paper* implies electron densities n e * 1 x 10 and B fields « 20 gauss 
for the case in which the X rays are primarily due to Compton interactions and 
the optical photons are due to synchrotron radiation* a rather higher value 
for n. and a lower value of B than Cavaliere and Morrison found for NGC 4151 
or 3C273. 

2. Thermal Bremsstrahlung and Blackbody Radiation 

The luminosity from an optically thin sphere emitting X-rays due to 

thermal bremsstrahlung at a temperature of T is 
-27 1/2 2 

L ' 2.4 x 10 T n. V erg/sec 

9 

where V is the volume of the source region and n fl is the electron density. 

For T * 10® K (as required by the spectral fits (Mushotsky et al. 1980') one 
requires n e * 1 x 10 14 era” 3 for V * 1 x 10 38 cm 3 . This gives an optical depth 
to Thompson scattering t » n at » 200. Therefore our initial assumption of 
optically thin emission cannot be correct and the output spectrum must be 
optically thick in shape if the X-rays originate in a sphere of hot gas. 
However the X-ray spectrum does not have this shape; therefore* we conclude 
that the X-ray emission without Compton amplification cannot be purely thermal 
in character. 


3. Therraa 1-Corap ton Models 
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Although simple accretion with thermal emieeion can explain the 
source of power# we need a more exotic emission mschanism to be consistent 
with both the spectrum and the rapid variability. Another X-ray source that 
shows rapid time variations in addition to a power law X-ray spectrum is Cyg 
X-l. We will now consider models originally constructed to explain Cyg X-1 
but scaled up to L ■ 10** ergs/sec. 

The more successful models have involved unsaturated 
Comptonization. For this mechanism to function the source must consist of a 
cloud of hot electrons with an optical depth to electron scattering less 
than a few and a copious source of soft "seed" photons. Under these 
conditions soft photons entering the cloud can be up-scattered in energy to 
X-rays# but do not spend enough time in the cloud to come into thermal 
equilibrium with the electrons. Shapiio et al. (1976) and Katz (1976) showed 
that the resultant spectrum is a power law for energies below kT e with an 
exponential fall-off at higher energies, tosdnyakov et al. (1979) and 
Takahara (1980) found qualitatively similar results in Monte Carlo 
calculations for a wide range of temperatures and optical depths. The 
location of the mil over is somewhat model dependent with the result from 
Shaprlo et al. occurring at the lowest temperatures. However# all of the 
calculated spectra roll over significally by 3 kT # . There is no evidence for 
either a roll over or a Wien peak in the spectrum of NGC 6814 or of any other 
AGN where spectra have been measured to higher energies. This implies that 
kT e > IS keV in NGC 6814. Since Can- A, the source with the beet determined 
spectrum# has a power law out to « 500 keV (Baity et al. 1981) it is possible 
that kT. > 300 keV. 

For unsaturated Compton models and kT* are related by 
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•• ' kT ' 


2/3 


(1) 


(Sunyaev and Titarchuk 1979/ , for T - 1.7 aa measured for HOC 6814. Using 
Equation 1 we compute that x • 4.0# 1.9# .8 for alactron taaparaturas of 30# 
100 and 300 kev. Notice that although kT # rangas over a factor of 10 tha 
product (t . kT # ) only changas by a factor of 2. This will bo useful below. 

In any mechai.ism involving inverse Compton scattering# tha harder 
photons are on the average scattered more times than ttu softer photons. 
Lightman at al. (19**8) pointed out that this implies that the harder photons 
spend more time in the source region# causing the more energetic photons to 
lag the soft during fluctuations. Lightman et al. suggested that the rise 
time of a flare t n , measured at energy E# should be given by 


t,, - t,* A ln(E/B,) (2) 

where t R ■ R/c (sise of source region), A is a constant on the order of unity 
and E ( is the energy of the original soft photons. More detailed calculations 
(Payne 1980 t L ghtman and Rybicki 1979) have shown that 

2 

ra c 

A " T kT *(3 “fiy (3) 

e 

where t ■ No^R and a is the power law index of the persistent 

source (a ■ r-1). Since (t . kT # ) appears as a product we can easily estimate 

A to within a factor of 2. If we use kT # « 100 kev# we find A « .6. 

In equation 2, t n is unobservable since it is measured from the 
(unobserved) injection time of the original soft photons. Me set ■ (l+f)t 0 
where t Q is the observed rise time of the flare and f is model dependent but 
only a slowly varying function of Eg. Solving for tg gives 
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U+f)t 0 

*R " .6 ln(B/B ) * 

8 


( 4 ) 


Before we can go on we mast estimate E Q . Lightman et al. (1978) 
suggested measuring the rise time of a flare in two energy windows and then 
determining the ratio of rise times 


D 



ln(E 1 /B a ) 

ln(B 2 /E s > 


(5) 


Solving for E s one finds that, for E s < .1 keV, E s is proportional to E^ and 
E 2 raised to large powers* Thus a small error (1%) in D can result in a 
factor of 10 uncertainty in E s . In addition, equation 5 is only approximate 
since it uses t m instead of the observable time t Q . An observational problem 
is that the statistics are not good enough to determine two different rise 
times in two narrow windows. Thus we are unable to determine E s accurately, 
but estimating this number will not result in too large an error in t R . 

Since Seyferts are known to be strong IR sources we choose E s *1 eV. 

This is consistent with the hypothesis of the soft seed photons originating in 
a large accretion disk which should be rather cool. He define t Q as the time 

for the intensity to go from .37 l MX to I . With log(E/E s ) * 3.8 we find 

f ' 2 using the function given by Payne (1980). Thus 

t ' 60 sec , 

R 

which says that the size of the source region is slightly smaller than the 
variation time. If we assume the peak flux comes from « 5 Schwarzschild radii 
then we estimate the mass of the compact object as 


' 9 ,F - 1 X 10* M . 

10 G O 


M - 
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If the geometry is canonical Karr then the peak X-ray flux cooes from 0.8 
Schwarz schi Id radii (Thorne 1974). For this case M * 8 x 10* Mq. 

Although it has been Impossible to measure D and thus check the 
consistency of the unsaturated Compton model , another test can be made. 

Photons at &2 (>Bj) should always lag photons at Bf by an amount 

At - t - t - At ln(E /E ) . (6) 

m m2 ml R 2 1 

This equation is independent of the soft photon energy but is a function of 
the electron temperature. It is possible to measure At^, using the 
cross-correlation function, for the entire data stream, not just during the 
rising portion of an event. Using tj ■ 60 sec and A •* .6 one ? .edicts At m * 

30 and 17 sec for the < 6 keV vs > 6 keV and MED vs HED 3 data respectively. 
Recall from the observations in Section XII that the measured times were -4 ± 
10 sec and 6 ± 3 sec, respectively. Both of these values of the lag time are 
consistent with the observed lag time being 1/3 or less of the predicted 
value. One possibility is that t R is smaller than our estimate of 60 sec (see 
below). Another point is that we are using rather broad energy windows which 
will "smear" the effect. Until we obtain more photons, allowing a smaller 
energy window, we cannot make any strong statements. However, the fact that 
the observed MED vs HED3 delay is in the right direction is encouraging. 

Both Payne (1980) and Lightman and Rybicki (1979) predict large (and so 
far unobserved) spectral changes during an intensity fluctuation. Me note a 
few reasons why we might not see such changes. First, Lightman and Rybicki 
suggested that one may identify the observed time for variability with the 
timescale for a change in the soft photon flux rather than tg. Equation (27) 
of Lightman and Rybicki is then applicable. This equation predicts n fi > 4 x 
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10 12 cm” 3 which would result In t R < 13 i«o« consistent with our uppar 
limit. However, if wa uaa tha datarmi nation of At ■ 6 i 3 aaoonda and apply 

A 

Equation 6 wa find t R « 22 i 11 aaoonda. Secondly, it ia poaaibla that tha 
output apactrum ca leu la tad by payna (1980) and Llghtman and Rybioki (1979) ia 
incorract bacaaua thay hava not calculatad a truly self -consistent modal. 
Therefore, wa cannot at praaant uaa our obaarvationa to critically conatrain 
such modala. Finally, it ia po88ible that unaaturatad Comptoniaation ia not 
tha corract modal. 

C. Source Lifetime 

We can aat an uppar limit for tha object' a lifetime at ita currant 
luminosity by requiring that it cannot hava accreted more material than ita 
praae.it mass at an assumed constant accretion rata. If we assume a 10% 
conversion efficiency of matter into radiation, consistent with soma of tha 
theoretical discussions of accretion onto a black hole (sea 'rises 1979 for 
discussion of the energetics of accretion onto black holes and sea Thorne 1974 
for a discussion of possible values of the efficiency), tha accretion rata for 
NGC 6814 must be > .002 M^j/yr to account for the observed luminosity. This 
gives a lifetime 


6 

10 M 


.002 M Q /yr 


5 x 10® yrs 


which is considerably shorter than the Nubble time. We cannot rule out the 

possibility that N > 5 x 10® Nq and that the efficiency is near 30% (Thorne 

1974), in which case the source could have maintained its present luminosity 

for a Hubble time. If the 10% efficiency aaaungition ia correct, one must 

consider alternate scenarios. Possibilities are that the object has a duty 

cycle of t /t . . . < .03 or that the high luminoeity is a recent 

L Hubble » 
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phenomenon. We note that « .02 of all galaxies with <-19 ara aotlva 
galaxies (Ruehra 1977). The valua of tha "on" duty oyola la tharafora 
oonaiatant with tha hypothaala that all galaxlaa show a Seyfert-like nature 
•or .02 of their lifetimes. 

D. Why Zs NGC 6014 Different? 

NGC 6814 la tha only AON In tha HKAO-1 A-2 sample (Tennant at al. 1981) 
to show such rapid variability. In this section wa speculate on why this may 
be so. For convenience wa have divided this section into two parts. First, 
we consider whether the same emission mechanism could account for the observed 
differences. Next we consider if we are observing a different emission 
mechanism but a "typical" galaxy. Throughout this section we will be using 
the fact that most active galaxies are ten times more luminous than NGC 
6814. However, if the observed variations are a property of luminosity only 
then both NGC 4151 and Cen A should show such rapid variability but do not, at 
least in the HEAO-1 data. 

1. Similar Emission Mechanism 

The light curve for NGC 6814 differs from other objects in two ways, 
the magnitude of the variations and the timescale involved. Let us now 
consider how decreasing the amplitude of variations would affect our ability 
to detect variations. For a jource counting rate of S and a background rate 
of B one can show 



k { (1 


£) -1 
s' ST 


1/2 


(7) 


2 2 2 

where T is the bin size, and k* • < o /a > is the average ratio of source 

6 C 

variance to the counting statistics variance* For the case where o is 

8 

2 2 

independent of time, the reduced x ie e measure of 1 + k. Thus equation 7 



1 ? 


2 

allows us to relate reduced x to the Intrinsic source variance. 

In order to minimise systematic errors we will compare NGC 6814 with 

other objects using only the counting rate from the small field of view of the 

xenon detector (which has only one- fourth of the total area). NGC 6814 data 

2 

with a 300 sec bin sise and total length of 2000 seconds gives x y " 10. Using 

B » 8.7 and S ■ 1.4 ct/sec we find that o /s » 36%. Other sources with 

s 

2 

counting rates similar to NGC 6814 typically have x y < 2 which implies 

that Og/8 < 1 3%. If other sources show similar variations then the amplitude 

must be reduced by at least a factor of three. A factor of three can be 

accounted for in a shot noise model. In this model one sees a superposition 

of identical shots each with height h and decay time T occuring at a rate A. 

2 2 

In this case S « Ahx and o « Ah t. For the case where h and T are constant 

s 

from source-to-source we find o /s = 1//A * 1//l. It will be possible to test 

8 

if this relation holds with our data for other sources. 

2 

A x test on a variable source will be indicative of variability if the 

timescale for variation is greater than the bin size and less than the total 

sample time. The HEAO-1 pointed observations are sensitive to all timescales 

from 10s of seconds to * 3 hrs. To observe longer timescales one needs to use 

scanning data. However, to obtain similar statistics a bin size of one to two 

days is required. Therefore timescales from 3 hr to 30 hr have not been 

adequately sampled, if we assume that the source variance is due to shot 

noise then it is possible that the brighter sources have the same shot rate 

but a longer decay time. In this case the more luminous sources would show 

longer timescales and again o /s * 1//L. 

s 

If rapid time variability only occurs in episodes, then our data 
indicate that at any one time only one object in « 30 will be having an 
episode. This can be tested only by watching a large number of objects in 
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addition to NGC 6C14 for a long period of time* 

2. Different Process 

In this section we assume that the X-ray emission from NGC 6814 is 
due to a different physical process than for other AGN's. We further assume 
that other sources also have a component similar to that in NGC 6814 but that 
for some reason it is generally not seen. Let us call the small variable 
component v and the larger more constant component C. We have to assume that 
the X-ray spectrum from v must be quite similar to the spectrum from C. 
Additionally we note that NGC 6814 is not unusual in its X-ray or optical 
properties. Thus the only obvious difference between V and C is the rapid 
time variability. 

Consider ways in which the V component could normally be hidden. 

One possibility is that the C component is normally much more luminous. If 

this case is correct then one might expect c /s a L . This is consistent 

s 

with the observation (Mushotzky et al. 1981) that on timescales > 12 hours the 
low luminosity sources are typically more variable. 

Another way that one might hide a V component is with geometry. This is 
possibly supported by the fact that NGC 6814 is almost a face on object. 
Heckman et al. (1978) surveyed 21 -cm emission from 23 Seyfert galaxies. Of 
this saoq>le no other active galaxy had a smaller HI line width than NGC 6814 , 
showing that most of the rotation is in the plane of the sky. If active 
galaxies send matter and or radiation out along the rotation axis then it is 
most likely that we would see it in NGC 6814. Even if we do not see such a 
"jet" directly then it might have blown a hole in the accreting gas 
surrounding the central object. This would allow us to see the nucleus 
directly whereas from other directions it would be obscurred by material 
accreting onto the central object. 
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VI. CONCLUSIONS 

NGC 6814 exhibits strong x-ray variability on timescales less than SOOO 
seconds. The autocorrelation function shows a characteristic time of * 100 
seconds. There does not appear to be a strong spectral-temporal effect since 
any change in the spectral index is less than .35 during strong intensity 
variations. There is a weak indication, however, that soft photons may 
systematically lead hard photons by - 10 seconds. 

He compare our data with various models and conclude that, unless 
reacceleration or continuous injection of particle occurs, neither the 
synchrotron nor the synchrotron self-Compton processes can be the origin of 
the X-ray emission. However, we do find t at our data are consistent with the 
general predictions of thermal-Compton models. He note that it is necessary 
to observe a roughly exponential rollover in the high energy end of the X-ray 
spectrum in order to confirm the validity of thermal Compton models of AGN 1 s. 

NGC 6814 appears to be unique. In our sample of over 30 active galactic 
nuclei (Tennant et al. 1981) only NGC 6814 shows rapid variability. He 
speculate that this might be related to its low X-ray luminosity. 
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APPENDIX A 

DATA CONSISTENCY CHECKS 

Since the rapid variability of NQC 6814 is so far unique we feel that a 
detailed discussion of the checks that we made will be useful. 

A. Instrumental Self Checks 

The Goddard Space Flight center's A- 2 experiment (Rothschild et al. 

1979) is a powerful tool to conduct stability studies of X-ray sources* The 
instrument consists of four detectors each with two different fields of 
view. Usually, three of the detectors (deck) view the source while the fourth 
(offset) detector views background. During the NGC 6814 pointed observation, 
each of the three deck detectors independently saw statistically the same 
light curve that is shown in Figure 1. while the offset detector had a 
counting rate consistent with a constant for the duration of the point. 

Instrumental stability is quite good and was tested during the so called 
"blank sky" points. A "blank sky" point occurred whenever the A-l instrument 
on HBAO-1. which looks 18 0 C from our look axis, was pointed at a source and 
there were no HBAO-1 sources (flux > 2 x 10“ 11 ergs/cm^ sec) in our field of 
view. Examining the look position from all A-1 pointings we found that 55 
different observations were truly "blank sky" for a total observation time of 
380 hours. We have analysed the data from these points and report that only 
once did we see a short duration rate increase in more than one detector. 

This event turned out to be time coincident with a solar flare. At no time 
did we see a short duration rate decrease. 

The rates during the "blank sky” points were analysed with our standard 
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data selection flags* At no tine did we see a clear indication of 
contamination- due to particle events. Therefore, we conclude that the net 
contribution due to particle contamination of the background subtracted 
counting rate is usually less than 0.1 count/sec in each detector. 

The "blank sky” points also showed that the background can differ by 4% 
from one point on the sky to another. Changing the background flux by 4% 
changes the average NGC 6814 flux by * 20% but does not significantly affect 
the shape of the light curve. 

B. Source ID 

Ariel-5 was first to report an X-ray source at the position of NGC 6814 
(Elvis et al. 1978). Their identification was uncertain since their error box 
also allowed the globular cluster PAL 11, which lies about 2° away, to be the 
X-ray source, rhe HEAO-1 error box (Mushotzky et al. 1980) was aligned in 
such a way that it is unlikely that the cluster is the X-ray source. For the 
pointed observation we can rule out the cluster with 100% certainty since it 
was never in our small (3° x 1.5°) field of view. For a picture of the two 
error boxes see Mushotzky et al. (1980). 

The Einstein Observatory observed the source with both the HRI and the 
I PC (Tananbaum 1980). The HRI showed a point source within 5" of the visible 
light nucleus which confirms that the galaxy is an X-ray source. The IPC 
showed no other sources with a flux greater than one- tenth the NGC 6814 flux 
in the 30' field of view (M. Elvis 1981). In addition, the IPC detected 
strong variability with a timescale of a few hours. Faster variations could 
not be detected since the count rate was only 0.07 ct/sec. We expect the IPC 
to detect a lower flux than HEAO-1 since it has a smaller area and it looks in 
a lower energy band which is affected by the low energy absorption seen in 
NGC 6814. We have converted the IPC counting rate into a spectral 
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normalisation by assuming a spectrum and using the tables in Oiacconi at al. 
(1979). Ms find that 0.07 IPC counts/sec correspond to e norm of 4.7, 2.2 and 
9.4 photons/cm* ksec keg for spectra with T ■ 1 .75, n h ■ 5.0} T ■ 1.50, Mg • 
2.9} r - 2.0, N h ■ 7.9 respectively where N H has units of 10 22 atoms/cm 2 . The 
three spectra represent the best fit, lower and upper 90% confidence range 
found by Hushotsky et al. (1980). The norm found by Mushotsky et al. for the 
total observation using the best fit values was 19 which is a factor 4 greater 
than the IPC norm assuming the same spectrum. However, if one compares the 
range of possible norms found for the IPC with the numbers in Tables III and 
IV we find that they do overlap. Thus the low IPC flux is not inconsistent 
with NGC 6814 being the X-ray source. 

C. Aspect 

HEAOl was designed as a scanning experiment and this has meant that the 
aspect solutions during a point are not as accurate as obtained during scan. 
The HEAO-1, A- 3 experimenters generate their own aspect solutions and have 
kindly provided some information to us. By comparing the two aspect solutions 
we find that the A-2 solution can drift up to 0.3° at the end of 6 hrs. Often 
the solution is much better. By looking at strong constant sources we have 
found that errors in the GSFC provided aspect solution can introduce an 
apparent 10% variation in source Intensity. He find that it is relatively easy 
to identify aspect problems since the area corrected rate is inversely 
proportional to the effective area. Thus although we can identify aspect 
errors we still cannot make any strong statements about source stability below 
the 10% level. 

For the NGC 6814 observation we have compared the observed flux with the 
effective area and see no correlation. He note that many of the variations 
are seen in the raw counting rate. Further we have conducted the stronger 
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teat of comparing tho counting rate with look position. Xf tha aapoot is in 
error or a confusing souros is inside tho detector beast one would aspect to 
see a simiiai flux whenever the spacecraft is looking at a similar positio in 
tho sky. Ws see no correlation between lode position and counting rate. 

During the rapid rate decreases in orbits me and three the spacecraft lode 
position was moving slowly through a region it had been looking at for more 
than 1/2 the orbit. 
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APPENDIX B 

N0H7SRX0DXC 8XGNAL ANALYSIS 

In this appendix m will present hm 4* tails of tho methods uitd to 
•nslyso the HOC 4814 light curve. No fool that such a proaoatation of this ia 
needed for oovoral reasons. First# thoro aro aovoral dofinitiona of tho 
autocorrelation function (ACT) and aovoral ways one can uao tho function. Wo 
diacuea our choice anf raaaona • Further# application of tho ACT to real data 
requires oonai do ration of the effocta of tho data gapa# finite data traina and 
photon counting atatiotica. Finally# it ia not trivial to coopute the 
uncertainties in tho derived results. This is due to tho fact that tho data 
are correlated and thus straight propagation of statistical errors can grossly 
underestimate the errors involved. For HOC 6814 we found that the dominant 
source of error arose from the finite observation length. Therefore we spend 
the most time discussing errors from this source. 

Given a binned set of observations Y^ with statistical errors we 
compute the ACF from 

p(u) - r(u)/r(0) 


where 


r(u) 


« < Y # * Y > < Y #^ * Y > 
. i N” u 1 l»u 

" N 1*1 o. o. 


i i+u 


B1 


where Y is the weighted mean and u is the lag time measured in bins. Jenkins 
and Watts (pl72f) show that using 1/N instead of l/(N-u) in the definition of 
r(u) introduces a bias into the estimator of the ACF for large u# but the 
biased estimator has lower mean square error. Equation si shows that the 
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2 2 

weighted ACT is directly related to reduced x • Thue x y should differ froe 

1.0 by eany sign* in order to determine a statistically significant ACT. This 

beoomss apparent as me goes to smaller bin sises since the statistical errors 

dominate and so becomes less significant. We have found that for practical 

2 

purposes the minimum bin sise occurs When x v differs from 1.0 by % 5 sigma 
«fh«rc sigma is given by /(2/M). For smaller bin sises noise dominates the 
ACT. 

Me now consider data gaps. Which can be of two types, short (a few 
minutes) or long ( » 1/2 hour). Short gaps have relatively little affect on 
the ACT, which means that any reasonable method of gap filling should work. Me 
gap fill by replacing the missing data with the average. This method will not 
introduce any false peaks into the ACT, since the ACT is only sensitive to 
deviations from the average. 

The long gaps present a somewhat different problem. One could again gap 
fill with the average but this is equivalent to assuming that the entire data 
stream is a sample of a stationary stochastic process. Me are uncertain 
whether this is a reasonable assumption. In particular, the large change in 
the average flux from orbit- to-orbit will increase the short '.erm 
correlations. If we assume that we are dealing with a stochastic process 
superimposed on top of a time varying mean, then tie can detrend the data to 
concentrate on the short term variations. Me "detrend” by subtracting the 
orbital mean from each orbit's worth of data. Another method to treat large 
gaps is to simply not cross the gap. This means that one computes the ACT for 
only one orbit at a time. Me prefer this method since it allows comparisons to 
be made between different observations. 

The dominant source of error in the timescales reported here is due to 


the finite length of the data stream. Crudely put, to measure a 



26 


timescale t on* n**ds an observation with a length of many (* 10) T. To 
eatimate the error Introduced by the finite data lengtn we havs generated 
e imula ted data streams using • fii:* 1 order autoregreeelve process (see Jenkins 
and Watts 1968). For a first order process a data point y(X) is generated by 


y(X) - oy(I-l) ♦ Z 
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where |a| < 1 and Z is a random number with mean <Z> and variance. var(Z). 
The mean and variance for the y's are given by 


<y> 

var(y) 


< 2 > 

1-0 

var(2) 
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The ACF for the first order process is an exponential with decay 
time t - -1.0/ln(a). Since the ACT is independent of the various moments of 
the data, the errors in the timescale will also be independent of the 
moments. The autoregressive process is useful since it is easy to compute and 
allows one to separate easily the timescale parameter a from the various 
moments. The shot noise model also allows this separation. 

The maximum likelihood estimate for a is p(1) (Jenkins and Watts 
1968). Although on* feels that it is possible to obtain a better result by 
using more than on* point, numerical tests have shown that the best estlnte 
is p(1) alone. It is possible to estimate the uncertainty in p(1) from 

o 2 » (1 *p(1) 2 )/L 
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(Box and Jenkins 1976) where L is the number of points in the longest# 
continuous data segment. He numerically tested the estimate for o for various 
L's and with a ■ .9 and believe the estimate is quite good for L > 100. To 
allow for photon counting noise# which is not considered above# and any other 
systematic effects the errors reported are 50% larger than the value given by 
B4. It is po *ble for p(1) +1 . 5o to exceed unity. This tells us that the 
timescale could be as long as the obsevation# which is reported as +* in Table 
II. 

For any finite observation of a varying signal# the ACF must show 
periods of negative correlations. He can estimate the lag time t z when the 
ACF crosses zero by finding the time when p(t) equals 1.5o. This gives 

t ' (T/2) ln(2N/3(l-p(1) 2 )) B5 

z 


For orbits one and two, N * 200 and T * 100 sec, equation B5 predicts t_ * 330 

z 

sec. The observed time for both orbits was 250 secs. This shows only that 
the zero crossing time is closely related to the observation length. 

Finally# Heisskopf et al. (1975) pointed out that photon counting noise 
affects the normalization of the ACF for non-zero lags. Their procedure for 
generating the correct autocorrelation function P c (u) was to multiply by a 
correction factor F(u) ; i.e.# p c (u) “ F(u) p(u). The function that we used 
has two values given by 


F(u) 



for 

u * 0 

for 

u * 0 


B6 


2 

where \ is the standard chi-square computed assuming the N data points have a 
constant mean. For NGC 6814# F was typically in the range 1.5 to 2.5. 
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Correcting (or photon noise means that the best estimate (or T is no 
longer p(l). To be consistent we determined t by making a least squares (it 
to ln(p(u) ) (or lags in the range 0 to t. This procedure is only slightly 
a((ected by data binning and (or the most* part gave estimates o( t consistent 
with p(l) . 

The spectral cross-correlation (unction CC(u) is based on equation Bl 
but now using two di((erent sets o( Y values. The cross-correlation lag time 
t is given by the ratio 

CC(1)/CC(-1) - exp(2t/T) (B7) 

if t is less than the bin size. Notice that the ratio CC(1)/CC(-1) does not 
need to be corrected (or photon noise and also is independent of a bias (to 
first order). To allow for the smaller expected error we apply equation B4 
with L replaced by N, the total number of data points used in computing the 
cross-correlation function. We use B4 to estimate the percentage error in 
p(l) and argue that the error in the ratio is twice the percentage error in 

4 

p(l). Since the ratio is independent of the bin size we used rather large 


bins to minimize Poisson errors 
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ORBIT 


1 


2 


3 


4 


TABLE 1 


NGC 6814 


MEAN FLUX 
cts/cm^ ksec 

5.6 

2.1 

4.8 

2.5 


2 

X 

(82 sec bins) 


279.6 


67.9 


77.2 


12.3 


DOF 


21 


23 


11 


9 


5 


2.9 


35.9 


6 
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ORBITS USED 


1 


2 


1-5 


1-5 


TABLE 2 


AUTOCORRELATION FUNCTION RESULTS 


DETREND 


TIMESCALE (SEC) 


NO 


85 


+60 

-25 


No 


230 


-130 


No 


248 


-110 


Yes 


101 


+59 

-26 


TABLE 3 




NQC 6814 

SPECTRAL FITS 


ORBIT 

A 

r 

H h (xlO 22 ) 

2 

X 

1 

20.4 

1 67 + ' 23 
1,67 -.21 

5 A +3 * 2 
5 ‘ 4 -2.2 

1.5 



+ .46 

. -+3.8 

4.8 

2 

7.2 


2 ‘ 5 -2.5 

3 

7.7 

+.23 

1 ‘ 35 -.21 

2 5 +2 * 5 
2 * -2.3 

14.0 

4 

5.7 

1 55 + * 51 
1,55 -.07 

0+3.6 

7.7 



+.52 

a +3.8 


S 

6.0 

• 

00 

1 

• 

o 

0 

3 • 
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TABLE 4 


HOC 6814 SPECTRAL FITS 


(ABSORPTION FIXED AT 4.3 X 10 22 ) 


ORBIT 


1 


2 


3 


4 


A 


16.5 


11.2 


11.1 


18.9 


1.59 + 


1.87 


1.49 


2.06 


2.07 


.07 

.07 


.19 

.16 


.10 

.09 


.24 

.21 


.27 

.22 


5 


20.4 
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FIGURE CAPTIONS 

Figure 1 - The X-ray flax from NGC 6814 ee e function of time for the first 
three spacecraft orbits* Data from both MED and HXD3 were added together to 
construct this curve, giving an effective bandwidth of 3-20 keV. The bin sise 
is 83 sec* Zero corresponds to April 28*00 (UT) , 1978* 

Figure 2 - a. Enlarged section of the NGC 6814 light curve plotted in Fig. 

1. Time of zero seconds corresponds to 24 minutes in Fig. 1. For this figure 
the bin sise was 20.48 sec. The data point at t * 500 sec occurred during a 
data drop out and so does not exist. 

b. The hardness ratios defined as > 6 keV flux divided by the < 

6 keV flux in HED 3 during the event shown in Fig. 2a. The hardness ratio was 
computed every 82 sec. The values of T for a given hardness ratio are 
indicated as straight lines labeled by the value of T. 

c. The counting rate for HED 1 during the event shown in Fig. 

2a. The HED 1 field was offset by 6° from NED and HED 3 and so was monitoring 

2 

the background. For the data shown x waa 42.06 for 42 degrees of freedom. 
Figure 3 - a. The ACF function for orbits 1-5. 

b. The natural logarithm of the ACF shown in 3a. 


c. The "detrended" ACF for orbits 1-5 
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d. The natural logarithm of 3a. 

Figure 4 • The photon index computed for each orbit ia plotted ae a function 
of the average flux during that orbit. We have fixed the absorption at Ng ■ 
4.3 x 10 22 atoms in the line of sight. The point near 9 x 10“'* erg/cn 2 
occurred during orbit 3 where the spectral fit was unacceptable. Thus the one 
sigma error shown may be underestimated. 

Figure 5 - The average hard flux is plotted whenever the soft flux was 
within a certain range. The line labeled r • 1.75 is the best fit line to the 
data. The fact that the line does not go through the origin indicates that 

the background subtracted was slightly in error, if the shape of the spectrum 

* 

is constant during an intensity change then the data should lie along a line. 
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